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The Renyi-Shannon entropy associated to critical quantum spins chain with central charge c = 1 
is shown to have a phase transition at some value ric ol the Renyi parameter n which depends on 
the Luttinger parameter (or compactification radius R). Using a new replica- free formulation, the 
entropy is expressed as a combination of single-sheet partition functions evaluated at n— depen- 
dent values of the stiffness. The transition occurs when a vertex operator becomes relevant at the 
boundary. Our numerical results (exact diagonalizations for the XXZ and Ji — J2 models) are in 
agreement with the analytical predictions: above Wc = the subleading and universal contribu- 

tion to the entropy is ln(L)(i?^ — l)/(4n — 4) for open chains, and ln(_R)/(l — n) for periodic ones 
(R=l at the free fermion point). The replica approach used in previous works fails to predict this 
transition and turns out to be correct only for n < nc. From the point of view of two-dimensional 
Rokhsar-Kivelson states, the transition reveals a rich structure in the entanglement spectra. 
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Introduction. — The entanglement entropy (EE) has 
become an important tool to probe and characterize 
many-body quantum states. In one-dimensional (Id) sys- 
tems, the celebrated logarithmic divergence of the 
von Neumann entropy of a long segment provides an efh- 
cient way to measure the central charge of a critical spin 
chain. In higher dimensions, EE can also be used to de- 
tect the presence of topological order [ ] , which is other- 
wise invisible to conventional local order parameters and 
correlation functions. Understanding how scales the EE 
of large subsystems has also opened a route to algorithms 
able to simulate efficiently these strongly interacting sys- 
tems in d > 1 [()]. 

In this letter, we use a seemingly completely differ- 
ent entropy to probe the ground-state of quantum spin 
chains, the Shannon entropy (or configuration entropy). 
For a normalized state lip) and a Renyi index n > it is 
defined as : 




where the states \i) form a basis of the Hilbert space. The 
basis states are chosen to be products of local states, and 
the Ising configurations = \ fti etc.) are the 
natural choice for the spin-^ systems we consider. This 
entropy was originally introduced as a mean to evaluate 
the EE in some particular two-dimensional (2d) states 
- the so called Rokhsar-Kivelson (RK) wave functions - 
and has allowed to investigate in details the EE at some 
some 2d conformal quantum critical points [ , ] . We will 
interpret our results in terms of 2d entanglement spectra 
at the end, but we will otherwise mostly concentrate on 
the Id spin chain point of view. 

A particularly interesting situation is that of the spin- ^ 
XXZ chain: 

■H = ^ (Sf ■ Sf_^_i + Sf ■ S'f^.i) + A ^ 5*^ • S'f+i (2) 

i i 



For a chain of length L, S'„ has a leading term propor- 
tional to L followed by universal subleading contribu- 
tions. If the chain is periodic, the first subleading term 
is of order 0(1) and was shown to be [7, 9] : 

5Pcriodic^(...)^^l^(^)_ Inn ^ (3) 

2(n — 1) 

where R is the compactification radius, related to the 
anisotropy A of the XXZ Hamiltonian (i?(A)^ = 2 — 
;| arccos(A)). With open boundary conditions the first 
subleading correction is also universal and takes the form 
of a logarithm of the length L of the chain [ ] : 

5,op™ = (...)i_lln(L), (4) 

which, in the 2d-RK language, confirms at c = 1 the 
arguments developed in [10] (from now on the term pro- 
portional to L will be omitted). 

In this letter we however show that these results are 
only correct below a critical value Uc of the Renyi pa- 
rameter. Using a replica-free formulation of the problem 
we prove that the Renyi parameter n effectively modifies 
the compactification radius of the chain (in a sense to be 
defined latter), and that a phase transition takes place 
at n = He when a (boundary) vertex operator becomes 
relevant. A central result concerns the location of this 
transition and the behavior of the entropy above n^. 

Uc ^ d^/R^ (5) 

Snynf'" = ^^-^ (nlni?-lno?) (6) 

«t =lnW^(f-i), (7) 

where d is the degeneracy of the Ising configuration with 
the highest probability Pmax in the ground state. At zero 
magnetization this configuration is c? = 2-fold degenerate: 

Nmax) = inn • • • n> or iw> = i;nt--nt>. (s) 
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Compact free field. — To understand this transition we 
adopt the 1+1 dimensional (Euchdian) point of view. At 
long distances, the model is described by a compactified 
"height field" h{x, r) with Gaussian probabilities: 



(9) 
(10) 



S[h] J dxdriyhy 

Z = j V[h] exp(-S'[/i]) 



where k is the stiffness, r the compactification radius 
[h = h + 27rr) and the physical Luttinger parameter 
(which fixes the decay exponents of the correlations func- 
tions) is i? = \f2Kr. For periodic (resp. open) chains 
h lives on an infinitely long cylinder (resp. strip) of 
perimeter (resp. width) L. In this language, the micro- 
scopic configurations \i) are replaced by configurations 
(\>{x) — h{x,T = 0) of the height field at r = 0. To eval- 
uate the probability p[(j)] we decompose the field h into 
an harmonic function h'^ which satisfies the boundary 
condition hf^^ = (j), and a "fluctuating" part Sh sat- 
isfying a Dirichlet boundary condition: h = h'^ + Sh 
[22]. Exploiting the Gaussian form of the action and 
Ah^ = 0, the "classical" and "fluctuating" part decou- 
ple: S[h] = S[h^] + S[Sh] and we get 



=exp(-5[/i*])- 



(11) 



where is the partition function of the whole cylinder 
(resp. strip) with a Dirichlet defect line at r = [23]. 
Now p[(/)] is raised to the (possibly non-integer) power n 



p[(j)]" = exp(-7iS'[/i'^]) 



z 



(12) 



and make the observation that exp(— nS" [h^]) is the 
Boltzmann weight in a system where the stiffness k has 
been replaced by k' = uk. From now on we explicitly 
keep track of the value of the stiffness k (as an index) 
and write: 



exp(-nS'K[/i'^]) 



cxp(. 



So, Eq. 12 can be written as 



(Pk[0])" =Pn« 



z. 



(13) 



(14) 



Inserting this result in Eq. 1 and using the fact that the 
probabilities p„k[(/)] are normalized, we get the main re- 
sult of this section: 
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nln 



(15) 



If we assume 2n = p to be an integer the partition func- 
tion Z{n) = X^-P? l"^^^ ^ natural interpretation in terms of 



p half-infinite systems glued together at their edge, form- 
ing a "book" with p sheets [] i]. The derivation above, 
however, never assumes 2n to be an integer and is there- 
fore different from the previous derivations involving a 
rephca trick [ , ]. Sn has been reduced to ratios of 
standard partition functions with respectively free and 
Dirichlet boundary conditions at the boundary between 
upper (t > 0) and lower (t < 0) parts of the cylindric (or 
strip). The complications associated to n-sheeted sur- 
faces [8, 9] have been avoided. For a compactified free 
field and the cylinder geometry, the ratio gjj = Zj^/Z^, 
is a well known "^-factor" [12, 13]: 



9l 



R- 



-1/2 



Combining Eq. 15 and 16 we find: 



InR 



Inn 



2(n-l)' 



(16) 



(17) 



in agreement with Refs. [7-9, 14]. For the strip geome- 
try, we also recover the n-independent result of Eq. 4 by 
applying the Cardy-Peschel formula [ ] to Eq. 15, with 
four angles 7 = 7r/2. 

Boundary phase transition. — Eq. 15 is a combination 
of two (7-factors and therefore probes the boundary of the 
system (the "bookbinding" in the book picture). The ac- 
tion at the lattice scale is not strictly Gaussian and other 
terms respecting the lattice symmetry and the periodicity 
h = h-\-2TTr are present, including vertex operators of the 
type cos(^/i). At the boundary such an operator renor- 
malizes to zero in the long distance limit if d^ > 2Kr^ 
[ l Otherwise it would lock the field to a flat configura- 
tion with degeneracy d. But the Eq. 15 cannot be valid 
anymore in this locked/massive phase since its deriva- 
tion assumes Gaussian probabilities. A boundary phase 
transition therefore takes place when the most relevant 
vertex operator become marginal in presence of a stiff- 
ness k' — UK [ ]. This gives the critical value of the 
Renyi index: 



d^ 
2Kr^' 



(18) 



In the case of the antiferromagnetic XXZ chain the two 
Ising configurations | titi ' ' ') and | |tit ' ' ') become 
the ground-states when A — >■ 00, which shows that an 
operator with two minima is present at the microscopic 
level. Taking 2Kr'^ = R{A) and d = 2, we find Eq. 5. We 
also note that the same argument applies to the Ji — J2 
(or "zig-zag") chain, and to any Luttinger liquid phase 
with umklapp terms ^ cos(2/i/r).[/' ] 

In this locked phase, the universal contribution to S„ 
is given by these d configurations only, so that 



1 



1-n 



In [d(p„,ax)"] 



(19) 



In the periodic case, Pmax simply corresponds to 
Z^ I Z and we recover Eq. 6. 
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h{0, r) = 7rr 



h{L^l.T) =7rr 



h{x,T = 0) = 7rr/2 

FIG. 1: (Color online) Height shift 5h — nr/2 for the semi- 
infinite strip with bottom boundary condition |imax)- 



The open chains turn out to be even more interesting. 
The universal contribution to Sn>na is encoded in 



Pn 



lim 



s e 



-rHl 



Z(U) 



(20) 



where -Z(U) is the partition function of a semi- infinite 
strip with bottom boundary condition limax) (see Fig. 1). 
The ratio in Eq. 20 is similar, but not identical to the one 
used below the transition (n < Uc). As before, four cor- 
ners with angle tt/2 will contribute to — Inpmax by a loga- 
rithm: — J InL. However, the configuration with highest 
probability does not exactly correspond to Dirichlet in 
the continuum limit since there is a height shift S be- 
tween the vertical edges of the strip, and the horizontal 
boundary (see Fig. 1). It can be treated by subtracting 
an harmonic function hs{x, r) = [25 /it) arg(x-l-ir), equal 
to on the horizontal boundary r = and 5 on the ver- 
tical boundary at a; = 0. The resulting contribution to 
the free energy is 



5F ' 



2^ 



InL. 



(21) 



The value shift 5 can be obtained using, for instance, 
a bosonization approach. The free boundary condition 
for the spins at the end of the open chain corresponds 
to Dirichlet for the free field and we have to set h{x = 
0,t) = h{x = L + 1, t) = irr to insure vanishing spin 
operators at both ends [17]. Then, the continuum limit of 
the configuration jzmax) corresponds to locking h{x, r — 
0) to the minima of the umklapp term cos(2/i/r) [ ], 
which has two degenerate minima in /i = 7rr/2 and 37rr/2. 
In both cases the height difference between the x — Q 
boundary and that at r = is (5 = 7rr/2. Summing 
up the contributions coming from Eq. 21 and from the 
Cardy-Peschel term, we recover our main result Eq. 7. 

Numerical simulations. — So far, for periodic chains, 
only the case n — 1 has been investigated numerically 
[7]. Fig. 2 shows the full n dependence of (constant 
term extracted by fitting the finite-size data) for different 
values of A. Although the system sizes are very small 
(L < 28) there is a good agreement with the theoreti- 
cal predictions, including the change of behavior at the 
predicted value of (which depends on A). It is only 
close to the Heisenberg point (A = 1) and above ric that 
the finite-size results deviate from Eq. 7. We attribute 




J2/Jl=0.2411 
A=0.4 
A=0.2 
A=0 

ln[R(A)]-log(n)/(2n-2) - 



FIG. 2: Constant term in the Shannon-Renyi entropy of 
periodic XXZ and Ji — J2 chains for different values of 
A and at the critical point J2/J1 = 0.2411. Each point 
comes from fitting the data for L — 20, 22, 24, 26 and 28 to 
aL + b + c/ L + d/ L'^ . Fat lines: theoretical prediction (Eq. 15). 
Eq. 3 is also plotted above ric (dashed lines) for comparison. 



these enhanced finite-size effects to the marginal opera- 
tors present at the SU{2) symmetric point. To circum- 
vent this difficulty we also studied the Ji — J2 chain at 
the critical value J2/J1 ~ 0.2411 [19], which has the same 
radius R = \/2 but where finite-size effects are much 
smaller. The data again agrees well with our prediction. 

Open chains were investigated in Ref. [ ] for various 
values of n and A, but the transition at ric was over- 
looked. In Fig. 3 we see a clear tendency for the loga- 
rithmic term to approach the theoretical curves, although 
the finite size data are still far from the thermodynamic 
limit. Interestingly, the entropies curves extracted from 
different system sizes cross in the immediate vicinity of 
the predicted value of ric- This lead us to conjecture 
that the coefficient of the logarithm may take a universal 
value at the transition point which only depends on the 
compactification radius R. At A = we pushed the nu- 
merics to L = 40 spins (see below) and the inset of Fig 3 
shows a collapse of the data for different system sizes 
onto a single curve in the vicinity of Uc = 4. This indi- 
cates a slow (^ L^^/^) but steady convergence towards 
a step function. We also get — —(1/6) InL with a 

high precision at A = 0, although we have no theoretical 
understanding of this value. 

A = m open chains. — This corresponds to free 
fcrmions and each probability pi is a determinant (Wick's 
theorem). Denoting by {xk} the positions of the up spins 
(fermions) and setting 9k = ■nxk/{L + 1) we have: 



v^ = p{{ek}) = 



L + 1 



L/4 



det (sin [j9k] 

l<j,k<n 



(22) 
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FIG. 3: Logarithmic term in the Shannon-Renyi entropy of 
open XXZ and Ji — J2 chains, extracted using a fit aL+b In L+ 
c + d/ L to four consecutive even systems sizes (L, • ■ • , L — 6) 
with L = 16 and L = 28. Fat lines: Eq. 15. Inset: scaling 
close to nc(A = 0) = 4. 

This determinant can be computed exactly [20]: 

t^'2sin2^-p. . 2 f 0j-Ok \ . 2 f 0j+(^k \ 
j=i fc>i ^ ^ ^ ^ 

and allows to go get Pmax: Pmax = 2^-^/^. The absence 
of a logarithmic term in S'n-j.oo ~ — Inpmax is consis- 
tent with Eq. 7, because the contribution from the height 
shift exactly compensates for the Cardy-Peschel term at 
i?(A = 0) = 1. This also rigorously confirms at A = 
the value of (5. We also checked the validity of Eq. 7 for 
the square lattice quantum dimer (2d-RK) wave function, 
for which d =\. In this case the microscopic height rep- 
resentation of dimer coverings allows to obtain b exactly 
at the lattice level. 

Conclusion. — The Renyi-Shannon entropy of critical 
Id systems does not only give a simple access to the Lut- 
tinger parameter of the system, it also shows that taking 
powers of the wave function gives rise to a whole line 
of critical points ending at a phase transition to an or- 
dered state. In fact, going back to the 2d RK point of 
view, this transition reveals a rich structure in the en- 
tanglement spectrum {Ei} of these 2d wave-functions. 
In some appropriate geometry [7] the probabilities pi are 
nothing else but the eigenvalues of the RK reduced den- 
sity matrix and the probabilities directly give the entan- 
glement spectrum Ei — — Inpj. The phase transition 
at ric shows that the spectrum has two sharply distinct 
regions in the thermodynamic limit: at high "energy" 
(small inverse "temperature" n) the universal contribu- 
tions to the probabilities are Gaussian whereas at low 
energy (large n) they are dominated by the ordered con- 
figurations. This contrasts with the fractional quantum 
Hall situation where the universal part of the entangle- 



ment spectrum only exist at low Ei [21]. These findings 
may open new directions in the study of quantum critical 
wave-functions - in Id and in higher dimensions - as well 
as their connection to boundary critical phenomena and 
f;-factors in particular. 
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